
NPS ARCHIVE 
1965 

GONCZ, J. 








v: 



A 



•.A. 



DUDLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOOL 
MONTEREY CA 93943-5101 



V 




V 





ULTRASONIC ATTENUATION IN 
SUPERCONDUCTING ZINC 



* * * i< i< 



Joseph P 0 Goncz 



ULTRASONIC ATTENUATION IN 
SUPERCONDUCTING ZINC 



by 

Joseph P. Goncz 

// 

Captain, United States Army 



Submitted in partial fulfillment of 
the requirements for the degree of 

MASTER OF SCIENCE 
IN 

PHYSICS 



United States Naval Postgraduate School 
Monterey, California 



19 6 5 



\r^> AV'c'A \)< 
1 -^ 0 ;": £ = f 




DUDLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOOl 
MONTEREY CA 93943-fmi 



ULTRASONIC ATTENUATION IN 
SUPERCONDUCTING ZINC 
by 

Joseph Po Goncz 

This work is accepted as fulfilling 
the thesis requirements for the degree of 

MASTER OF SCIENCE 
IN 

PHYSICS 
from the 

United States Naval Postgraduate School 



ACKNOWLEDGEMENTS 



I wish to thank Professor John R 0 Neighbours for his enthusiastic 

guidance and patient assistance in conducting the experiment, I also 

wish to express my appreciation to Mr, Lynwood F, May, USNPGS, who built 
3 

the He cryostat, for his assistance and also to Lt. David Dickmann, USN, 
who wrote the computer programs. 



ii 



ABSTRACT 



Measurements of the ultrasonic attenuation of 10 and 30 Mcs/sec 

longitudinal waves by pulsed-echo techniques were made on a 99.999% 

pure single crystal of superconducting zinc in the |^000lj direction 

as a function of temperature from 4.2 K to 0.320 K using an open-ended 
3 

type He cryostat » The attenuation was found to be frequency dependent 
and decreased less sharply near the superconducting transition tem- 
perature, T , than predicted by the Bardeen-Cooper-Schrief f er (BCS) 
theory c Attenuation due to electron-phonon interactions only was 
found by subtracting from experimental points the value of residual 
attenuation gotten by extrapolation of the data to T = 0°K. Using 
the BCS theory the zero degree superconducting energy gap was found 
to be 2 e(0) = (3.36 + 0„13)kT c with T c = 0.817°K. 
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INTRODUCTION 



The rapid decrease in the attenuation of ultrasonic waves in 

superconductors below the superconducting transition temperature, T , 

12.. C 

was first observed by Bbmmel and Mackinnon . 5 Bommel investigated 

the ultrasonic attenuation of shear and longitudinal waves in a cylin- 
drical crystal of lead in the frequency range from 9-27 Mcs/sec from 
approximately 18°K to 1.5°K. He found a marked decrease in the atten- 
uation in the superconducting state which was frequency dependent. No 
theoretical explanation was made, though he suggested that the results 
might lend support to the idea that superconductivity is due mainly to 
interactions between electrons and lattice vibrations. Mackinnon in- 
vestigated the absorption of 10 Mcs/sec sound pulses in a superconducting 
polycrystalline tin rod and in a polycrystalline magnesium rod between 
1.5°K and 4.2°K. In tin he found the attenuation decreased by a factor 
of about 50 and in magnesium that the attenuation decreased steadily 
with decreasing temperature. 

This decrease in the electron contribution to the ultrasonic atten- 
uation was first satisfactorily explained in 1957 by the theory of 

3 

Bardeen, Cooper, and Schrieffer (BCS theory). The BCS theory predicts 

that the ratio of the electron contribution to the attenuation in the 

superconducting state, O' , to that in the normal state, O' , is 

s n 



O' 

s 

a 

n 




where 2e (T) is the temperature dependent energy gap. The theory further 
predicts that as T approaches zero, the energy gap approaches 3.50 kT^ . 

4 

Morse and Bohm were the first to apply this theory to the analysis 
of ultrasonic attenuation data on tin and indium, shortly before the 
publication of the BCS theory, basing their analysis on a conversation 
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with Professor Bardeen. Their analysis pointed out the remarkable 

agreement between the BCS theory and their experimental data. Since 

5 6 

then the BCS theory has been applied to many superconducing metals. * 

These investigations have all shown close agreement between the 
theory and experiment, but have also indicated the need for refinement 
of some of the simplifying assumptions of the BCS theory. The theory 
shows that the basic interaction causing superconductivity is that of 
the exchange of a virtual phonon between two electrons. An electron 
moving through the lattice creates a virtual phonon which is absorbed 
by a nearby electron. If the distortion of the lattice is such as to 
surround the electron by a positive charge greater than the electronic 
one, the second electron sees an attractive potential greater than the 
repulsive Coulomb one and is attracted. These electrons form what is 
known as Cooper pairs. The theory then shows that the lowest energy of 
the system is achieved by the maximum pairing of electrons of opposite 
spin such that their momenta add to zero. The difference in the normal 
ground state energy and the superconducting ground state energy is due 
entirely to an interaction potential, less than zero, arising from the 
correlation between the Cooper pairs. The energy of this interaction 
leading to the scattering of one pair to another is the difference 
between the energy of interaction by means of a phonon and the energy 
of Coulomb interaction. 

The existence of an energy gap then follows from the small, but 

finite, amount of energy required to break up a Cooper pair. But more 

important, having a single unpaired electron removes a large number 

of pairs which could have occupied the state of a Cooper pair. This 

gap will be temperature dependent since increasing numbers of pairs 

will be broken up by thermal excitation as the temperature rises from 

0°K to T . 

c 

Morse, Olsen, and Gavenda^ pointed out the need to consider the 
anisotropy of the energy gap and showed that in tin there is considerable 



- 2 - 



anisotropy (3.2 kT^ - 4.32 kT^) . In order to conserve both momentum and 
energy in a phonon-electron interaction, the component of the electronic 
velocity parallel to the direction of sound must equal the phonon veloc- 
ity; i.e., the sound velocity. Since the velocity of electrons at the 
Fermi surface is in general much greater than that of sound, only those 
electrons moving at essentially right angles to the direction of sound 
propagation can interact with the phonons. On the Fermi surface these 
electrons lie in a narrow band perpendicular to the sound velocity. 

Since the ultrasonic attenuation in a particular direction is an average 
over this narrow band, the energy gap measured will be an average for 
the electrons lying within the band. Anisotropies in the Fermi surface 
will then lead to anisotropies in the energy gap. The BCS theory makes 
the simplifying assumption that the interaction energy is isotropic and 
constant near the Fermi surface. It is this assumption that leads to 
the BCS prediction that the zero degree energy gap is the same for all 
superconductors and is equal to 3.5 kT^ . 

The relation for the ratio of ultrasonic attenuation in the super- 
conducting state to that in the normal state follows from the coherence 
of paired electrons. In the normal state scattering of an electron to 
another state is independent of any other transitions. However, the 
pairing of electrons in the superconducting state makes these transitions 
interdependent and the theory shows that the scattering by pairs is 
destructive. Thus, the absorption of phonons, causing ultrasonic atten- 
uation^ decreases sharply below T^. With the assumption of low energy 

phonons; i.e., hw « 2e(0) , the ratio & /o' follows. 

s n 

The present work on zinc, which has a hexagonal crystal structure, 
was undertaken to determine the energy gap in the direction of the 
c~axis and to compare the results with measurements made in the base 
plane. Ultrasonic frequencies on the order of 10^ cps were used, 
well satisfying the condition for the given expression of Q^/o^ * 
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EXPERIMENTAL PROCEDURES 



A cylindrical, single Crystal of zinc, approximately 2 \" long and 
3/4 n in diameter with the c-axis, [OOOl] direction, parallel to the 
cylinder axis was purchased from Semi-Elements, Inc . The purity of the 
sample was stated to be in excess of 99.9997<>. From this ingot the experi- 
mental sample was prepared to the approximate desired size, limited by 
the cryostat, by cleaving both ends and etching. The ingot was held in 
a styrofoam form, immersed in liquid nitrogen and cleaved by laying a 
new, sharp, razor blade at the desired point and rapping it with a 
mallet. To reduce the diameter, the ends were coated with brushing 
lacquer and the ingot etched in dilute HC1 acid. 

Standard ultrasonic pulse-echo techniques^ were used and will be 
discussed only in terms of the experiment. The success of this method 
depends strongly on having the end faces of the sample flat, parallel, 
and perpendicular to the direction in the crystal one wishes to investi- 
gate. To achieve this, the zinc sample was cast with a cerro-alloy in 
a bored out, cylindrical jig whose faces were machined flat, parallel 
and perpendicular to the crystal axis to a tolerance of 0.003 n . The 
ends were first lapped on metallographic paper and then polished on a 
diamond dust wheel. 

The final sample was (3.1602 + 0.0002) cm long and approximately 
5/8" in diameter at room temperature. A Laue back reflection pattern 
of the zinc crystal showed concentric circles, indicating a thin poly- 
crystalline layer on each end. However, the x-ray spots showing through 
the circles indicated that the faces were perpendicular to the c-axis, 

A V, X-cut, 10 megacycle quartz transducer was next bonded to the 
zinc sample. The success of this bond is a critical step in the experi- 
mental procedure. After some experimenting, the following procedure 
proved to be reliable. Both the quartz and sample faces are cleaned first 
with acetone to remove greasy substances, second with alcohol to remove 
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water, and finally with acetone to ensure the surfaces are dry. The 

best bonds were obtained on days when the humidity was low. Next a 

very small amount of analytical grade glycerin was applied to the zinc 

and the transducer centered and squeezed to make the glycerin flow and 

also to remove air bubbles. Glycerin proved to be the most satisfactory 

bonding material, though its use required the making of a new bond for 

each run. The critical temperature for the bond was from 40°K to 50°K. 

If the bond held on being cooled through this range; the bond was good 

for the run. However, as can be seen, one has already started to cool 
4 

with He and if the run must be stopped for a broken bond, the delay is 
both costly and time consuming. 

A schematic for the ultrasonic pulsing equipment is shown in Fig. 1. 
Radio frequency pulses (10 Mcs/sec and 30 Mcs/sec) of approximately two 
microseconds' duration were generated by an Arenberg Pulsed Oscillator 
triggered by a Time Mark Generator with a repetition rate of 1000 cps 
or 100 cps. The slower rate had to be used below 0.5°K to reduce heating 
of the sample. The rf electrical pulses were simultaneously fed to a 
Type RM 45A Tektronix Oscilloscope, externally triggered by the Time 
Mark Generator, and impressed across the piezoelectric quartz transducer. 
The zinc furnished the rf ground and the oscillator, tuned to the funda- 
mental or odd harmonic of the transducer, caused it to vibrate and send 
sound waves through the zinc. The sound waves reflected off the other 
end return to the transducer to produce a voltage across the electrodes 
furnishing the initial pulse and also to be reflected back through the 
sample. With a pulse rate of 1000 cps, these reflected echoes had 
adequate time to decay away between pulses. The voltage produced by each 
returning echo was amplified and displayed in proper time relation to 
the initial pulse on the oscilloscope. The change in height of each echo 
is a measure of the attenuation in the sample due to phonon-electron 
interaction, impurities, imperfections and in the transducer bonding 
system. Absolute attenuations are gotten by taking pictures of the scope 
traces with a Polaroid Oscilloscope Camera, Model 196A, Special, using 
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ULTRASONIC PULSING EQUIPMENT SCHEMATIC 
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3000 speed/Type 47A Polaroid film and plotting echo heights vs echo 
number on semi- log paper. Attenuation effects other than those due to 
phonon-electron interactions are subtracted off as discussed in the 
treatment of the data. 



Data were taken in the temperature range from 4.2°K to approximately 

o 3 

0.320 K, The temperature range was obtained in a He cryostat of the 

open-ended type, a schematic of which is shown in Fig. 2. It is desirable 

to cool the sample as fast as possible to reduce cooling time, but this 

is tempered by the necessity of not straining the bond. 



The sample was loaded in a copper sample holder, inserted in the 
3 3 

He cryostat and the He chamber flushed with dry nitrogen and left at 

3 

one atmosphere of dry nitrogen The vacuum jackets of the He cryostat 
4 

and the He dewar were thoroughly flushed with dry nitrogen and left 

4 

soft at approximately 300 microns dry nitrogen pressure. The He dewar 
was also flushed and left at one atmosphere dry nitrogen pressure. The 
nitrogen dewar was then filled with liquid nitrogen and the whole system 
cooled to approximately 77°K which takes about six hours. 



The He and He jackets were then pumped out until the pressure was 

- 6 3 

less than 5 x 10 mm Hg . Keeping the pressure less than this in the He 

jacket was essential to prevent thermal conduction at the lowest temper- 

4 

atures. While this was being done, the nitrogen and He traps were 
3 

filled and the He gas circulated for approximately one-half hour to 

4 

condense out all impurities except He . This was found to be necessary 

as impurities would freeze out in the cryostat and around the sample 

3 4 

holder causing throttling when pumping full out on the He . The He 

impurity caused no trouble as thermometry was done using a carbon resistor 

3 

and not by using the He vapor pressure. 



When the jackets were pumped out and the He purified, the He dewar 

3 

was filled and about eight mm Hg of He let into the cryostat; the re- 
mainder being returned to the reservoir. The system was then cooled to 
4.2°K, taking about 12 hours. The process could be accelerated by 
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HELIUM 3 REFRIGERATOR SCHEMATIC 




Fig. 1 



8 




softening the He jacket with He However, achieving a satisfactory 
vacuum after the system was cooled down was very difficult since helium 
at that temperature is very difficult to pump. Without a good vacuum 
the low temperatures desired could not be reached. 

o 4 3 

Once 4,2 K was reached, the He was pumped on and the He let 

into the cryostat. A minimum temperature of 1.09°K can be reached 

3 

this way, well below the temperature at which He at one atmosphere 

will condense; i.e., 3ol95°K. When equilibrium was reached, the cir- 

3 

culating pump was valved off to pump on the He . At 500 micons the 

3 

booster pump was cut in and pumping full out on the He the lowest 
temperatures were reached. 

Temperatures were determined by measuring the resistance of a 

Vwatt Speer Carbon Resistor, nominal resistance 470 ohms, 107o, grade 

1002, in a Wheatstone Bridge Circuit employing a lock-in amplifier to 

measure the bridge unbalance. The resistor was calibrated against the 

4 o 4 

vapor pressure of He down to 1.1 K using the M 1958 He Scale of Temper- 
ature". Down to the lowest temperature reached, the resistor was cali- 
brated against the magnetic susceptibility of cerium magnesium double 
nitrate (CMN), Ce^Mg^ (NO^) ^ *24H^0 using a Cryotonics Mutual Inductance 
Bridge. ^ The susceptibility of CMN is strongly anisotropic with the 

susceptibility parallel to the symmetry axis essentially zero. To obviate 
the necessity of packing single crystals, powdered CMN was used. The 
effective susceptibility is then approximately 65% of the value perpen- 
dicular to the symmetry axis, but the gain in the packing fraction off- 
sets this decrease. 

Since the temperature dependent part of the susceptibility of CMN 
perpendicular to the symmetry axis follows Curie's Law 

%■? 

the mutual inductance bridge readings are correlated to the temperature 
and thereby calibrating the carbon resistor. 
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EXPERIMENTAL RESULTS 



Measurements of the ultrasonic attenuation of 10 and 30 Mcs/sec 
longitudinal waves were made in the [0001] direction on a single crystal 
of zinc whose ends were slightly polycrystalline . Data were taken in 
the temperature range 4,2°K to approximately 0»320°K and in the earth's 
magnetic field. The superconducting transition temperature was determined 
to be T^ = 0.817°K and the zero degree superconducting energy gap was 
measured as (3.36 + 0.13)kT^* 



To measure the attenuation as a function of temperature, pictures 
of echo trains were taken at different temperatures and the heights of 
the echoes plotted on semi- log paper to obtain the attenuation in units 
of attenuation per echo. These were plotted as a function of temperature 
and when extrapolated to T = 0°K gave an estimate of the residual atten- 
uation. Subtracting the residual attenuation from experimental points 
gave the attenuation due solely to electron- phonon interactions. These 
data were then analyzed using the BCS function for reduced attenuation 
vs reduced temperature 



a* = 



exp 



sill 

kT 



+ 1 



The details of all these analyses are given in Appendices I - III. 
In order to illustrate the methods, Figs. 3-14 graphically display the 
results for one set of data; those of 13 November 1964. 



After the transition temperature was determined and an estimate of 
the residual attenuation made, computer program RESATT was used to 
determine the residual attenuation and the zero degree energy gap. 

Figs. 3-5 show the point plots of reduced attenuation vs reduced 
temperature for three slightly different values of residual attenuation 
subtracted from the experimental data. By inverting the BCS expression 
the following function Y is defined: 

Y = l/ln(2/0f* - 1) = (kT /6(t))t. 

c 



- 10 - 




X -SCALE - 2.00E-01 UNITS/INCH. 
Y-SCALE - Q.60E-01 UNITS^IHCH. 



GQNCZ PROJECT 8321 13 NOU 64 38MCS 
REDUCED ATTENUATION US REDUCED TEMPERATURE 



- 11 - 





m 002 00H 006 009 010 

X-3CALE - Q.00E-01 UHIT3/IHCH. 

-SCALE - Q.00E-0J UNJTS/'JNCH. 

GONCZ PROJECT 0321 13 NOU 64 30MCS 
REDUCED ATTENUATION US REDUCED TEMPERATURE 



Fig. 4 



- 12 - 





X -SCALE - 9.88E-8J miTS/'JNCH. 

Y -SCALE - Q.88E-81 miTS'INCH. 

GONCZ PROJECT 0321 13 NOU 64 30MCS 
REDUCED ATTENUATION US REDUCED TEMPERATURE 



Fig. 5 



- 13 - 




X “SCALE - Q.ME-ei LIMITS'' INCH. 

Y “SCALE - Q.00E-ei UNITS/ JMCH. 

CONCZ PROJECT 0321 13 NOU 64 30MCS 
Y-l/(LOG(Q/ALPHA-l» US REDUCED TEMPERATURE 

. ~ Fig. 6 



- 14 - 




X -SCALE - 2.0SE-0J UNJT3/JMCH. 

'f-SCALE « Q.08E-QJ UNJTS/'JNCH. 

GOliCZ PROJECT 8321 13 NOU 6 1 ! 30MCS 
Y=l/(L0G(2/flLPHfi*-l)) US REDUCED TEMPERATURE 

_ Fig. 7 



- 15 - 




I 




X-3CALE - Q.ftdE-01 UNJTS^JNCH 
'f-JCWLE - O.0QE-81 UNJT3/JNCH. 

CONCZ PROJECT 0321 13 NOU 64 30MCS 
Y= l/(L0G(2/flLPHfl * ~ D) US REDUCED TEMPERATURE 

Fig. 8 



- 16 - 





X “SCALE - Q.08E-01 UNJTS/'JNCH. 
Y-SCALE - Q.eeE“8 1 UNITS^JMCH. 



GONCZ PROJECT 8321 13 NOU 61 30MCS 
REDUCED ATTENUATION US REDUCED TEMPERATURE 



Fig. 9 



- 17 - 





aaa m em 006 000 010 



K “SCALE - Q.80E-0J UNITS/ INCH 
Y-SCALE • 2.00E-0J UNJTS/JNCH. 

GONCZ PROJECT 0321 13 NOU 64 30MCS 
REDUCED ATTENUATION US REDUCED TEMPERATURE 



Fig. 10 



- 18 - 





K-3MLE * Q.00E-01 LMJT3/JNCH. 

Y-3CflLE - Q.00E-01 UNIT SCINCH. 

GONCZ PROJECT 0321 13 NOU 64 30MCS 
REDUCED ATTENUATION US REDUCED TEMPERATURE 



Fig. 11 





X-3MLE - Q.68E-01 UNJTS'JNCU 
'f-3CflLE - Q.00E-01 LJMIT3/JMCH. 

GONCZ PROJECT 0321 13 NOU 61 30MCS 
Y=l/(L0C(2/flLPHfl*-l)) US REDUCED TEMPERATURE 



Fig. 12 



-20- 



I 





X -SCALE - Q.08E-01 UNJT3/JHCH. 

'f -SCALE - Q.ME-ei UNITS/'JHCH. 

CONCZ PROJECT 0321 13 NOU 61 30MCS 
Y=l/(L0G(2/RLPhfl*-l)) US REDUCED TEMPERATURE 

_ _ _ __ Fig. 13 



- 21 = 





X -SCALE - 2.06E-01 UNIT5/JMCH. 

'f -scale - Q.eeE-ei units/jnch. 

GONCZ PROJECT 0321 13 MOU 61 30MCS 
J'CfsVfl_0G(2/ALPHfl * ' D) US REDUCED TEMPERATURE 

Fig. 14 



- 22 - 




Figs. 6-8 depict the point plots of the function Y vs the reduced 
temperature, t , corresponding to the same values of residual atten- 
uation. All these figures have the corresponding curve plots for the 
BCS expression for 2e (0) = 3.50 kT^ . 

The criterion used to find the residual attenuation and the zero 
degree gap was that on the plot of Y vs t the curve best fitting the 
experimental points should become a straight line passing through the 
origin for t less then 0.35. This follows from the expression for Y 
which becomes linear for 2 /o'* much greater than 1.00. The results of 
these analyses are given in Table I. 

Once the residual attenuation was determined, computer program 
VARIGAP was used to show how changing the zero degree energy gap would 
change the BCS curves for the residual attenuation and Y, Figs. 9-11 
show the BCS curve for reduced attenuation vs reduced temperature for 
a given residual attenuation and with changes of 0.10 kT^ in the energy 
gap determined from RESATT. Figs. 12-14 give the corresponding results 
for Y vs t . As can be seen the BCS curves do not shift much for the 
value of incremental change used. 

It was found that the attenuation was constant above T^ and frequency 
dependent. Table I shows that attenuation increased with increasing 
frequency. The higher attenuation for 30 Mcs/sec in the normal state 
on 5 January 1965 compared to that of 13 November 1964 reflects the 
greater attenuation losses from the poorer bond. The value for the 
zero degree gap was gotten from computations based on semi- log plots. 

It is felt that this method is better than that of comparing echo 
heights since the latter method involves propagating through all the 
data any error made in measuring the echo height chosen as the basis for 
comparison. However, Table I shows that essentially the same results 
are obtained using either method. The difference in the methods is 
apparent in the different values of residual attenuation obtained. Taking 
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the value at the lowest temperature reached for comparing echo heights 
required that the residual attenuation be determined from this value 
rather than from an average of low temperature values as was done 
using semi-log plots. 
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TABLE 1 



RESULTS OF RESATT COMPUTATIONS 



Attenuation from Semi- log Plots 





Run 




Frequency 

Mcs/sec 


Normal 

Attenuation 

a db/cm 
n 


Residual 
Attenuation 
a Q ,db/ cm 


2£±0X 

KT 

c 


Results 

a , db/cm 
o’ KT 

c 












0.310 


3,46 






13 


Nov 


64 


30 


3.092 


0.320 


3.58 


0.320 


3.58 












0,330 


3.62 
















0.525 


2.96 






5 


Jan 


65 


30 


3.350 


0.565 


3.26 


0.565 


3.26 












0,605 


3.49 
















0,660 


2.73 






5 


Jan 


65 


10 


1.376 


0,700 


3.25 


0.700 


3.25 












0.740 


3.63 













Attenuation from Comparison of Echo 


Heights 








Run 




Frequency 


Normal 


Residual 


2£(0) 


Results 








Mcs/sec 


Attenuation 

a db/cm 
n 


Attenuation 

a , db / cm 
o 


KT 

c 


a , db/cm 
o 


2C(0) 

KT 

c 












0.290 


3.45 






13 


Nov 


64 


30 


3.092 


0.300 


3.58 


0.300 


3.58 












0.310 


3.65 
















0.525 


2.75 






5 


Jan 


65 


30 


3.350 


0.565 


2.99 


0.585 


3.30 












0,605 


3.57 
















0.700 


3.01 






5 


Jan 


65 


10 


1.376 


0.715 

0.730 


3.33 


0.715 


3.33 
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DISCUSSION OF RESULTS 



The value obtained for the transition temperature, T^ = 0.817°K, 
is lower than those of other investigators, However, those other 
experiments were determinations of the zero magnetic field transition 
temperature for zinc and show considerable variation,, In heat capacity 
measurements, Phillips^determined T^ = 0,825°K, Geballe and Matthias^ 
in work on isotope effect measured T = 0,855°K and Cochran and 

ll ^ Q 

Mapother in critical field experiments found T^ = 0,875 K. If a 
parabolic relationship between reduced magnetic field and reduced 
temperature is assumed 



H 

H 



1 




y 



where is the critical magnetic field above which zinc will not 
become superconducting, then the transition temperature is reduced 
about five millidegrees for = 50 gauss and H = 0,5 gauss. Thus 
T^ = 0o817°K agrees favorably with the lowest value quoted above. 



Table II gives the results of several investigators who have 
measured the zero degree energy gap in zinc. There have been no ultra- 
sonic attenuation measurements along the c-axis other than this one. 

The result of this work is higher than the value found by Biondi et al 
in a microwave absorption experiment in which there were "violent 

oscillations" in the data for hv/kT less than 2,7, The measurement 

c 

by Zemon and Boorse can only be taken as an average over all directions 
in the crystal. Bohm and Horowitz found the attenuation so great 
along the c-axis that they were unable to make any measurements. How- 
ever, their measurement of the energy gap in two directions in the base 
plane and the measurements of the energy gap along the c-axis by 
Biondi and this work show anisotropies in the energy gap. These results 
also indicate that the gap along the c-axis is less than that predicted 
by the BCS theory and also is less than that for directions in the 
base plane. 
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TABLE II 



ZERO DEGREE ENERGY GAP IN SUPERCONDUCTING ZINC 



2 (0)/kT 


Direction 


T 

c 


Method 


Investigator 


3.8 + 0.2 


[l21o] 


0.85 


Ultrasonic 
Attenuation 
233 Mcs/sec 


Bohm and 
Horowitz 


3.4 + 0.2 


[ io To] 


0.85 


Ultrasonic 
Attenuation 
233 Mcs/sec 


Bohm and^ 
Horowitz 


3.00 + 0.15 


[oooi] 


0.838 


Microwave 
Absorption 
15-100 Gcs/sec 


Biondj 
et al 


3.01 + 0.09 


Polvcrystalline 

Sample 


0.825 


Microwave 
Absorption 
50-75 Gcs/sec 


Zemon ^jid 
Boorse 


3.36 + 0.13 


[oooi] 


0.817 


Ultrasonic 
Attenuation 
10 & 30 Mcs/sec 


Present 

Work 
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A more precise determination could be made if lower temperatures 
could be reached. The lowest temperature reached in this work was 

0. 320°K; i.e., 0.38 t. Thus to determine the energy gap, the experi- 
mental data had to be extrapolated some 40% on the reduced temperature 
scale . 

The graphs of reduced attenuation vs reduced temperature can be 
fitted with a BCS curve quite well for low t. However, for t approaching 
1.00, the experimental points lie above the BCS curve. Bohm and Horowitz 
don f t report the comparison of their points with the BCS curve. How- 
ever, investigators of other elements have found just the opposite; 

1. e., points for higher t lie below the BCS curve. For Niobium, 

Weber ^ found the effect to be frequency dependent. With T^ = 9.15°K, 
the BCS curve for the measured value of the energy gap, 3.63 kT^ , 

fit the experimental points for 220 Mcs/sec longitudinal wave along 

the [l00] axis. However, for 30 Mcs/sec the points lay below the BCS 

curve. He also found an isotropic result for the energy gap in the 

three crystallographic directions and suggested that any anisotropy in 

the gap may have been smeared out by impurity scattering. For thallium, 

T^ = 2.38°K, Saunders and Lawson^ found that attenuation of longitudinal 

waves decreased more sharply below T^ than predicted by the BCS theory 

and that the anisotropy in the energy gap ranged from (3.76 - 4.10)kT . 

c 

These and similar results on other elements are in contrast to the 
findings of this work which should be investigated further. One basic 
assumption tacit in all these computations was that the residual atten- 
uation was temperature independent. Though measurement above the 
transition temperature bore this out, the attenuation in the normal 
state for temperatures less than the zero magnetic field T^ should be 
measured. This might be done by placing the sample in a magnetic 
field to destroy the superconducting state. 

Finally, measurements at higher frequencies should be carried out. 

Bobetic^ has carried out the calculations giving the ration of 0t lot 

s n 
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according to the BCS theory for an isotropic model at arbitrary 
frequencies and temperatures and predicted a discontinuity in attenu- 
ation at frequencies corresponding to the energy gap. This would 

2 

require frequencies of the order 10 Gcs/sec, 

In this work attempts were made to make measurements at 50 and 
70 Mcs/seCo However, so much power was needed to get an echo, that 
the sample heated excessively. Using a shorter sample would help 
solve this problem. 
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APPENDIX I 



DATA ANALYSIS 

Analysis begins by measuring the echo heights on each picture. 

This is most accurately done on an X-ray film reader to + 0.01 mm. These 
heights are then plotted on semi-log graph paper vs echo number; one 
graph per picture. The slope of this graph is absolute attenuation per 
echo. 

The basic problem in determining the absolute attenuation in the 
temperature range attained (4 0 2°K to approximately 0.320°K) is that the 
attenuation is so great from slightly below the superconducting transi- 
tion temperature, T> to 4 0 2°K that one usually has only one echo. 

c 

Below T c , the attenuation drops rapidly and enough echoes appear to 
obtain the attenuation graphically. 

A typical plot of absolute attenuation, a, vs temperature, T, is 
as follows: 
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where 



a = absolute attenuation 

Oi =: residual attenuation 

0 

Oi = attenuation due to conduction electrons in the superconducting 
state 

(X = attenuation due to conduction electrons in the normal state 
N 

a. - attenuation determined from raw data 

1 

T = Temperature - °K 

T = Superconducting transition temperature 
c 

A Oi = Change in attenuation between the superconducting and normal states 



The 
db/cm is 



slope of the semi-log plots are in echo 
given as follows: 



The attenuation in 



Attenuation (db/cm) 



Slope of Semi-log Plot x 8.6858 
2 x sample length in cm 



We assume an exponential decay of echo heights 0 



A = A 1 e 



-Oi* n 



Let 

th 

A^ - amplitude of n echo in the superconducting state 

A = amplitude of n ^ echo in the normal state 
N 

A 1 = amplitude of the initial pulse 
n = echo number 

OiT = attenuation in the normal state from T to 4.2°K; a constant 
c c 

as verified experimentally 

Oi - temperature dependent attenuation in the superconducting state 
L = sample length in cm 



Then 



A n = A'* 



• n 



A. = A e 

l 



“Q^-n 



Ln A v , = Ai A' - n a_ 
N T 



Ln A. ss Ai A' - n CC. 

l l 



in(-i)= n (a - a i ) = - nAa 
N c 
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or 



GL 



/,,/ , 8.6858 0 S'i , 

(db/an) = 5^ U <—) + a t 

N 



To determine OL from the single echo pictures above T , one com- 
T c c 

putes with n = 1, for as many point by the above equation as one 

T c 

can from the semi-log plots* The best horizontal straight line fit to 



a plot of these points gives OL 



Tc' 



Then using the equation 



0^ (db/ cm) 



8 c 6858 
2nL 




+ a v 



gives the attenuation for the remaining points below T c where there 
are not enough echoes to determine the attenuation from the semi-log 
plots . 

Extrapolating the plot of OL. and CL vs T to T = 0 will give an 

1 *c 

estimate of OC . Where the curved part intersects the value for O' 

o T c 

determines T • 

c , 0L n 

ie S 

The quantity desired is the reduced attenuation OL = — — , the 

N 



ratio of attenuation due to conduction electrons alone in the super- 
conducting and normal states vs t = T/T^, the reduced temperature. 
Subtracting from OC ^ and OC ^ the quantity OC^ which is attributed to 
impurity scattering*, end losses, transducer losses, etc., gives 



* 



OC 




oc. - a 

1 o 




Subsequent calculations are very sensitive to the value 0^, which 
is not easily determined since one must extrapolate approximately 407 o 
on the reduced temperature scale to t = 0 for zinc. 

The BCS theory gives 

* 2 
a = ^€(T)/kT + l 

for the variation of OC for the longitudinal attenuation with T, where 
€ (T) is the temperature dependent energy gap between normal and super- 



conducting electrons and k is the Boltzmann constant. What we want is 
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e (0) o The above expression can be inverted to give 




This is the equation of a straight line passing through the origin 
•k 

when 2/a » 1. Hence, the slope at t = 0 for a plot of y vs t gives 

the reciprocal of the zero degree superconducting energy gap in units 

of kT o 
c 

Data can also be analyzed by comparing the echo height of the 
same echo on each of the pictures. Then for a given n 

_ , „ , x 8c 6858 o , A no x 

“i (db/c " > = 2 SE — U ( T ~- 1 + “o • 

ni 

One must have some reference attenuation here denoted a^ to signify 
the attenuation at the lowest temperature reached. The attenuation can 
be most accurately determined at this point since the maximum number of 
echoes is present. 
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APPENDIX II 



PROGRAM RES ATT 

The raw data; i,e M pictures* must be first analyzed in the manner 
outlined in Appendix I, Then the computer can be utilized to compute 
the reduced attenuation and the function Y, used to determine the zero 
degree superconducting energy gap* and plot these values separately as 
a function of the reduced temperature. Program RESATT* written in 
Fortran 60 language and compatible with the 1604 CDC computer, will 
accomplish these calculation. The primary purpose of this program is 
to determine the residual attenuation and to find the zero degree energy 
gap. The residual attenuation will be that value which when subtracted 
from the experimental values of attenuation will make Y, for t less than 
0.4, plot as a straight line intercepting the origin. 

For each set of data the program provides for the point plotting 
of computations based on three chosen values of the reduced attenuation, 
each set plotted on a separate graph. On each of these plots a theo- 
retical curve based on the BCS expression, using a chosen zero degree 
energy gap, can be graphed. This is computed from values of £(T)/6(o) 

vs t, the ratio of the temperature dependent energy gap to the zero 

18 

degree energy gap versus reduced temperature as follows : 

* 2 

a = C( DEGAP 
2t 

/ + 1 

Presently 44 values of C(I) are used. More can be computed from the 
reference • 

The title* which must be changed for each set of data, identifies 
the experimenter, date on which data were taken, frequency, and computer 
project number. These are within the program itself. The program notes 
follow. 

Code 

TTR(I) - Reduced temperature, t, for theoretical BCS plot 
ATS (I) - Reduced attenuation, CC , for theoretical BCS plot 
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"k - 1 

TY(I) - (ln(2/a -1)) , for theoretical BCS plot 

M - Number of data points for theoretical BCS plot, maximum number, 100 

EGAP - Zero degree energy gap, 2 £ (0) , (no kT c ) 

C (I) - 6(TTR(I)/6 (0) 

N - Number of experimental values of T & 0^ (up to 60) 

TC - Experimental value of Transition Temperature 

AN - Experimental value of a = a 

n T« 



AZ(0) - Three experimental values of OL ^ 

T(I) - N experimental values of temperature 
A(I) - N experimental values of attenuation 
RT(I) - N values of reduced temperature 



AS(1, I) 
AS (2, I) ( 
AS (3, 1) ^ 



N values of 

* a - a 0 

a = — 

a - a 

n o 



for ( AZ(0) - small 

^AZ(O) - medium 
(^AZ(O) - large 



AY(1 , 1)^ 
AY(2, 1) 1 
AY(3,I) 

/ 



N values of 
1 




for 



AZ(0) 
AZ (0) 
AZ(0) 



- small 

- medium 

- large 



Plots: Two sets of point plots 

k 

First: Three plots of OL vs t, 

Experimental Points AS(1,I) vs RT(I) for AZ(0) small - □ 

Experimental Points AS(2,I) vs RT(I) for AZ(0) medium - O 

Experimental Points AS(3,I) vs RT(I) for AX(0) large - 

Each will have the theoretical BCS curve ATS(I) vs TTR(I) plotted 
on the same graph. 

Second: Three plots of Y vs t, 

Experimental Points AY(1,I) vs RT(I) for AZ(0) small - D 

Experimental Points AY(2,I) vs RT(I) for AZ(0) medium - 

Experimental Points AY(3,I) vs RT(I) for AZ(0) large - 

Each of these will have the theoretical BCS curve TY(I) vs TTR(I) 
plotted on the same graph. 
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For each run a total of 3 + M + N cards are needed in the following 
order. The last digit of each imput in each column in the stated order 
is to be in the 10, 20, 30, etc., column of the data card. 



Number of cards 
1 

1 - M 
1 

1 - N 
1 



Imput 

M, EGAP 

Card No., TTR(I) , C(I) 

N, AN, TC 

Card No., T(I), A(I) 
AZ(1), AZ(2), AZ(3) 



Following is a complete print out of RESATT with the card deck 
composed for the data of 13 Nov 64. Card 000-197 comprise the program 
and cards 200-243 are the MUhlschlegel values. As outlined, the cards 
to be changed for a different set of data would be 57-61, 199, and 
244-265. 
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APPENDIX III 



PROGRAM VARIGAP 

Once the best value of the residual attenuation and a determination 
of the zero degree energy gap has been made from RESATT, VARIGAP will 
compute and plot the theoretical BCS curves of reduced attenuation vs. 
reduced temperature and the function Y vs reduced temperature for three 
different zero degree energy gapso On each of these will be the cor- 
responding experimental points for one set of data and the chosen re- 
sidual attenuation. To the extent that they differ from RESATT, the 
program notes for VARIGAP fol.low 0 



CODE 








EGAP (I) 


- Three values of the zero degree 


energy gap (no kTc) 


AS 


- One value of the residual attenuation 




ATS (1,1) 
ATS (2, I) 
ATS (3, I) 


) M theoretical values of 

7 reduced attenuation from 

J the BCS expression 


for * 


^EGAP - small 
\EGAP - medium 
EGAP - large 


TY(1 , 1)') 
TY(2,I)> 
TY(3, I ) j 


M theoretical values of 
Y from the BCS expression 


for | 


^EGAP - small 
J EGAP - medium 
EGAP - large 


AS (J) 


- N experimental values of reduced 


attenuation 


AY (J) 


- N experimental values of Y 






PLOTS : 


Two sets of curve plots 







First: Theoretical reduced attenuation vs reduced temperature 

ATS (1,1) vs TTR(I) for EGAP - small - □ 

ATS (2, 1) vs TTR(I) for EGAP - medium -<} 

ATS (3, I) vs TTR(I) for EGAP - large - A 

Each of these will have point plotted on them the N experimental 
values of AS(J) vs T(J) . 

Second: Theoretical Function Y vs Reduced Temperature 

TY( 1 , I) vs TTR(I) for EGAP small 
TY(2,I) vs TTR(I) for EGAP medium 
TY(3,I) vs TTR(I) for EGAP large 
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Each of these will have point plotted on them the N experimental 
values AY(J) vs T(J). 

For each run a total of 3 + M + N cards are needed plus the correct 



title cards. The 


last digit of each imput in each column in the stated 


order is to be in 


the 10, 20, 30, etc., column of the data card. 


Number of Cards 


Imput 


1 


M, EGAP(low), EGAP(medium) , EGAP(large) 


1 - M 


Card No., TTR(I) , C(I) 


1 


N, AN, TC 


1 - N 


Card No., T(J), A(J) 


1 


AZ 



Following is a complete print out of VARIGAP with the card deck 
composed for the data of 13 Nov 64„ Cards 000 - 188 comprise the program 
As outlined, the cards to be changed for a different set of data would 
be 57 - 61, 190, and 235 - 257. 
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APPENDIX IV 
PROGRAM ZEROGAP 

This program provides for a least square fit of experimental points 
to the BCS function 



a 

s 




2 



C(T) 




1 



by minimizing the square of the vertical distances from the curve to 
the experimental points • This is done using two parameters; the re- 
sidual attenuation and the zero degree energy gap. The output of the 
program is the value of these two parameters » Since in the BCS ex- 
pression 0^ and a involve only the attenuation due to electron-phonon 
interaction, all other causes of attenuation must be subtracted out. 
Thus, 



a* -a 



a -a 
N o 



-fin 

„ k T ! 
1 + 1 



or 



a' = a + (a XT -a ) 

o v N o 



Hi i 



& 



kT 



+ 1 



where ex' and a„ would be on a BCS curve shifted an amount a due to 
N o 

residual attenuation. 

What is desired is to minimize the residue 

N 

R = Y. (a-a ') 2 

i=l 



where a is an experimental measurement of attenuation and N is the number 
of experimental points. 

Let the following quantities be defined which will be used in the 
program p 

OL - A(I) - Experimental measurement of attenuation 

t - RT(I)- The reduced temperature corresponding to the temperature T(I) 
at which A(I) was measured. 
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OL - AZ - Residual attenuation 
o 



one of two parameters to be determined 
C(I)-US(I) - The Miihlschlegel values of £(T)/ 6(0) vs t 0 



By use of an interpolating program, the values of C(I) were obtained 
from the Miihlschlegel values for 0.191 4 t 4. 1.000 in increments of 
0o001 and incorporated in the program. In the printout of the program 
following they are cards 200 - 910, where for example, 200 means 
t=0o200 and the number in the 10th column is C(I) for the t. 

C(I)/t - UC(I) 

The equation for CL ' can then be written 



CL 1 = CL + (CL -CL ) — — ■■■ ■ ■■- 

o + ' N o J C(I)° 2 6(0) 

^ 2t 



+ 1 



Now let FER(I) = 



C(I) • 2 6(0) 

2t EXP 

£ +1 



UC(I) -EGAP 
2 



+ 1 



where EGAP is the zero degree energy gap = 2 6(0) - the other parameter 
to be determined. Then the residue becomes 



R = 



N 



I 

i=l 



A(I) -AZ- (AN-AZ) 2FER(I) 



2 



From this an explicit expression for AZ 
condition 



<3R 

<2AZ 



0. 



can be gotten from the 



Carrying out the steps gives 



AZ 



[ (2ANFER(I)-A(I) j |^2FER(I)-l] 
^ [ 2FER(I) -1 j [2FER(I)-l] 



in which terms cannot be canceled since they are under the summation 
signs. In this expression only FER depends on EGAP 0 
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The other condition for minimizing the residue is 



<3R 

£k EGAR 



0 . 



Carrying out the steps gives 



3R 

d EGAP = 



N 



2 I 

i=l 



jA(I)-AZ-2(AN-A2)FER(I)j 



X 



[2(AZ-AN) FER(FER-l) ] 



So what is needed is to minimize £ |^(2ANFER(I) -A(I) ) -AZ(2FER- 1) J UC 

X £ FER(FER-l)j j^AN-Azj 

subject to the above value of AZ, The remainder of the program notes 
are the following : 

ONE (I) = 2FER(I) - 1 

TWO(I) = 2ANFER(I) - A(I) 

The imputs needed are the best guess of EGAP; AN, the experimental 
determination of the attenuation in the normal state; TC, the transition 
temperature; and the N data points of A(I) vs T(I) . 

One card contains N, AN, TC, and EGAP, the last digit of each is in 
column 10, 20, 30 and 40 respectively. Next follow the N data cards of 
A(I) vs T(I) ; one card for each set with the last digit of the card number 
in column 10, that of T(I) in 20, and that of A(I) in 30. A print out of 
the program follows at the end of the appendix „ The imput cards follow 
the blank card which must be added to the end of the program cards. 

Table III gives the results of program Zerogap and compares them 
to those of program Resatt. Figs. 15-26 give the plots of the BCS curve, 
using the values of residual attenuation and the energy gap determined 
from Zerogap and the point plots of the corresponding experimental data. 

It is evident that the program does what it is designed to do. 
However, the values of residual attenuation are markedly different in the 
30 mcs/sec data from that obtained by extrapolation of the curve of 
attenuation vs temperature. For the 30 mcs/sec, 5 Jan 65, semi-log data 
the value obtained is unrealistic in that it is greater than most of the 
data points. Only the attainment of lower values of reduced temperature, 
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t, can resolve the question as to whether or not the value of residual 
attenuation obtained from Zerogap is realistic* 

The program seems to work best when there are data points over 
the whole range of t; e.g., 30 Nov 64 and 5 Jan 65, 10 mcs/sec data. 

Also, when the points are not too scattered as in the 30 Nov 64 data, 
the results are reasonable 0 The values of (2*85 and 2.82)kT c obtained 
for the energy gap for these data are reasonable bearing in mind the 
previous observation that for higher t the experimental points lie above 
the 3.50 kT^ BCS curve. As one can see the lower the value of the energy 
gap, the higher the BCS curve. 

This program can be of use in interpreting data if one realizes 
that the analysis depends completely on the BCS theory and the simplifying 
assumptions used in obtaining the expression for the reduced attenuation 
and masks the possibility of other phenomenon taking place. 
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TABLE III 



COMPARISON OF ZEROGAP AND RESATT RESULTS 



RUN 




FREQUENCY 


TYPE 


ZEROGAP 


RESATT 








Mcs/sec 




RESIDUAL 


ENERGY 


RESIDUAL 


ENERGY 












ATTENUATION 


GAP 


ATTENUATION 


GAP 












db/cra 


2€(0)/kT c 


db/cm 


26(0) /kT 


13 


Nov 


64 


30 


Semi-log 


0.1772 


2.8544 


0.320 


3.58 


13 


Nov 


64 


30 


Echo Heights 


0.1507 


2.8171 


0.300 


3.58 


5 


Jan 


65 


30 


Setni-log 


0.8755 


7.2968 


0.565 


3.26 


5 


Jan 


65 


30 


Echo Heights 


0.4387 


2.5304 


0.585 


3.30 


5 


Jan 


65 


10 


Semi-log 


0.7166 


3.5941 


0.700 


3.25 


5 


Jan 


65 


10 


Echo Heights 


0.6097 


2.1181 


0.715 


3.33 
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X -SCALE - 2.Q8E-01 UNIT.VINCH. 

Y-SCALE - 2.60E-Q1 UHITS/IMCH. 

GONCZ ■ PROJECT 0321 13 NOU 64: 30MCS 
REDUCED ATTENUATION US REDUCED TEMPERATURE 



Fig. 15 
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> -SCALE - L.C0E-01 UNITVINCH. 

Y -SCALE ** 2.OOE-01 UNITS' INCH. 

GOMCZ PROJECT 0321 13 NOU 64 30MCS 
Y-ML0G(2/ALPHA*-1)) US REDUCED TEMPERATURE 



Fig. 16 
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S -SCALE - -aGOE-OI UNITS' INCH. 

Y -SCALE - 2X0E-01 IMTS^IMCH. 

GOMCZ PROJECT 0321 13 NOU 6^ 30MCS 
REDUCED ATTENUATION US REDUCED TEMPERATURE 



Fig. 17 
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X -SCALE - Q.80E-01 UNITS/INCH. 

Y -SCALE - O.08E-O1 L1MITS/IACH. 

GOMCZ PROJECT 0321 13 NOU 64 30MCS 
Y=1/(L0G(2/ALPHA*-1)) US REDUCED TEMPERATURE 



Fig. 18 
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X “SCALE - 2.8CE-Q1 UNITS/INCH. 



Y “SCALE - Q.8SE-01 UNITS/INCH. 



GONCZ PROJECT 0321 5 JAN 65 30MCS 
REDUCED ATTENUATION US REDUCED TEMPERATURE 



Fig. 19 
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X-SCfiLE = 2.001-01 UMITS^INCH. 
Y-SCfiLE = 2.00E-01 UNIT^IMCH. 



GONCZ PROJECT 0321 5 JAM 65 30MCS 

Y= ML0G(2/ALPHA * - D) US REDUCED TEMPERATURE 



Fig. 20 
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x -scale - 2.005-01 units, inch. 

Y -SCALE = 2 .CSE -01 UNITS. INCH. 



GO.NCZ PROJECT 0321 
REDUCED ATTENUATION US 



5 JAN 65 30MCS 
REDUCED TEMPERATURE 



Fig. 21 
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. vvru. i_ i j t ji 

Y-SCPLE - S.CSE-01 UNITS- INCH. 

GOMCZ PROJECT 0321 5 JAN 65 30MCS 

Y-l-'(L0G(2/RLPHfl*-l)) US REDUCED TEMPERATURE 



Fig. 22 
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'-SCALE - tk.QCE-01 UNITS' ’INCH. 
’-SCALE - '2.Q2E-01 UNITS^INCH. 



CONCZ PROJECT 0321 
REDUCED ATTENUATION US 



5 JAN 65 IONICS 
REDUCED TEMPERATURE 



Fig. 23 
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X-SCPi-E - 2J32-01 UNITS'-INCH. 

Y-ecaE - 2.ecz-0i umitvimch. 



GOMCZ PROJECT 0321 5 JAN 
Y=--l.'CL0G(‘2/flLPH,q*-l)) US REDUCED 



65 18MCS 
TEMPERATURE 



Fig. 24 



68 - 




X-535LE - teeZ-BI LMJ75VMH 
Y-SJfiLE - MBE-01 IMTS'IMCH. 



G0MC2 PROJECT 032 1 
REDUCED ATTENUATION US 



5 JAM 65 10MCS 
REDUCED TEMPERATURE 



Fig. 25 
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: -KfiLE - 2.KE-01 UtITMNCH. 
Y -SCALE - 2.OSE-01 U'ilTS/IMCH. 



GOMCZ PROJECT 0321 5 JAN 
Y-ML0G(2/flLPHR*-l)) US REDUCED 



65 10MCS 
TEMPERATURE 



Fig. 26 
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9 CONTINUE 00000310 

PRINT lO*(T(I)*A(I)»RT(I)tUC(I)»I=l»N) 00000320 

10 FORMAT ( 4E20 * 8 ) - 00000330 

PX=EGAP ' ' 00000340 

TRY=EGAP-0.1 00000350 
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APPENDIX V 
PROGRAM AVERGAP 



This program provides for a least squares fit to the BCS function 
using one parameter; the zero degree energy gap with the residual atten- 
uation specif iedo Thus the program and the program notes are the same as 
those of Program Zero-gap, except that the residual attenuation is part 
of the imput data rather than being derived from an analytic expression . 

The program will compute the zero degree energy gap which minimizes the 
square of the vertical distance from the experimental points to the BCS 
curve for three different values of residual attenuation * Thus s the 
expression 

£ £ (2ANFER(I) - A(I)) - AZ(2FER(I) -l)j UC [ FER(FER-lj J^AN-Azj 
i=l 



is minimized for each given AZ. 

The type of data cards, their format, and order of imput are the same 
as for Program Zerogap, with the exception that one card need be added. 

This is the last card and contains the three values of residual atten- 
uation; the last digit of each to end in column 10, 20, and 30 respectively. 
These values are to be determined as in Program Resat t; i.e., from an ex- 
trapolation to T = 0 on a plot of the experimental data of attenuation vs 
temperature . 

Table IV gives the results of Program Avergap and also a summary of 
all the computer calculations . It is seen that the results of Avergap 
are more consistent with those of Resatt than are those of Zerogap and 
again reflect the fact that the experimental points for higher t lie above 
the BCS curves used in Resatt . The results for the data of 5 Jan 65, 30 
Mcs/sec still seem unreliable with such a high energy gap. This is due to 
the fact that there are data points for low t only. 
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RESATT AVERGAP ZEROGAP 

RUN FREQUENCY RESIDUAL ENERGY GAP RESIDUAL ENERGY GAP RESIDUAL ENERGY GAP 
Mca/sec ATTENUATION 2€(0)/kT c ATTENUATION 2£(0)/kT c ATTENUATION 2£(0)/kT c 
db/cm db/cm db/cm 
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TABLE IV. SUMMARY OF COMPUTER CALCULATIONS 
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